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Abstract 

The good performance of the GW approximation for band-structure calcu- 
lations in solids was long taken as a sign that the sum of self-energy diagrams 
is converged and that all omitted terms are small. However, with modern com- 
putational resources it has now become possible to evaluate self-consistency 
and vertex corrections explicitly, and the numerical results show that they 
are, in general, not individually negligible. In this review the available data is 
examined, and the implications for practical calculations and the theoretical 
foundation of the GW approximation are discussed. 

1 Introduction 

Many-body perturbation theory jlj represents a powerful method for studying the 
properties of interacting electron systems from first principles that goes beyond the 
limits of local mean-field approaches. It is based on the Green function, which 
can be interpreted as a propagator that describes the evolution of an additional 
electron or hole injected into the system and interacting with its environment via 
the Coulomb potential. In this way the Green function can be directly related to 
experimental photoemission spectra in solid state physics, and macroscopic observ- 
ables like the total energy are obtained through an integral with the respective 
quantum-mechanical operator. Furthermore, unlike the Kohn-Sham eigenvalues in 
density- functional theory |2[ |j| , the resonances of the Green function have a well 
defined physical meaning as electron addition or removal energies and correspond 
to the proper quasiparticle band structure. They can be calculated from a modified 
single-particle Schrodinger equation, in which exchange and correlation effects are 
rigorously described by the so-called self-energy operator, which takes the form of 
a nonlocal, energy-dependent potential Q. 

A practical way of calculating the self-energy is given by a perturbation expan- 
sion in terms of the Coulomb potential and the Green function of the corresponding 
noninteracting system jjj. This expansion is most conveniently written in the lan- 
guage of Feynman diagrams §, which may be related to distinct scattering mecha- 
nisms. As it is not possible to sum the complete infinite perturbation series, physical 
intuition can serve to identify the most important contributions. Of course, the se- 
lection of diagrams depends on the nature of the problem. In simple metals and 
semiconductors correlation is predominantly long-range, because electrons interact 
with their environment through polarization of the surrounding medium and thus 
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avoid proximity. This effect is well described by Hedin's GW approximation H, 
which includes dynamic polarization in the random-phase approximation. 

The performance of the GW approximation has been discussed in several re- 
cent reviews [0,^]. Most importantly, it corrects the systematic underestimation 
of semiconductor band gaps in Kohn-Sham density-functional theory, giving very 
good agreement with experimental data Another more subtle effect relates 

to the narrowing of the occupied band width in the alkali metals, which is also 
described accurately Because of this success it was long believed that the 

GW approximation indeed captures all relevant self-energy terms for these materials 
and that the excluded diagrams have negligible weight. However, recent numerical 
studies that explicitly evaluated self-consistency and vertex corrections, the prin- 
cipal omissions in the GW approximation, have cast doubts on this assumption 
and suggest that the high quantitative accuracy instead stems from a cancellation 
of errors. A better understanding of these effects would not only strengthen the 
theoretical foundation of the GW approximation, but might also show the way to 
further systematic improvements. For this reason I review recent progress in this 
field below and assess the available data. 

This paper is organized as follows. In Sec. || the GW approximation is in- 
troduced. Self-consistency and vertex corrections are discussed in Sees. || and ||, 
respectively. Finally, the main points are summarized in Sec. ^[ Atomic units are 
used throughout. 



2 The GW approximation 

Many-body perturbation theory is based on the Green function, which is defined as 
the expectation value of the time-ordered operator product jlj 

G(l,2) = -i(¥|T$(l)^(2)]|*>. (1) 

The short-hand notation (1) = (ri,<7i,ti) indicates a set of spatial, spin and tem- 
poral coordinates, \^) denotes the normalized ground-state wavefunction in the 
Heisenberg picture, and T is Wick's time-ordering operator that rearranges the 
subsequent symbols in ascending order from right to left. Besides, ■ijj '(2) and ^>(1) 
represent the electron creation and annihilation operator in the Heisenberg picture, 
respectively. In the absence of a time-dependent external potential, the Green func- 
tion only depends on the difference t\ — t% and can be mapped to frequency space 
through a one-dimensional Fourier transform. In the case of noninteracting systems, 
where the wavefunction is a single Slater determinant, this yields the expression 

GoCn.r*;*,) = £ ^f. 1 ^"^^ , (2) 

^ u - e«k + «sgn(e„ k - 

in terms of the solutions 

-^V 2 + V s (r)j <^„ k (r) = e„ k (^ nk (r) (3) 

of the single-particle Schrodinger equation. Here /i denotes the chemical potential 
that separates occupied from unoccupied states and r\ is a positive infinitesimal. 
Spin variables have been suppressed, because in the absence of magnetization Gq is 
symmetric and diagonal in a . In the following it is assumed that the single-particle 
potential V s already includes the Hartree potential. 

The Green function of an interacting electron system is related to the propagator 
of the corresponding Hartree system through Dyson's equation |l3| 

G(l, 2) = G (l, 2) + [g {1, 3)E(3, 4)G(4, 2) d(34) . (4) 
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The self-energy operator X rigorously describes all exchange and correlation ef- 
fects. It can be expanded in a perturbation series comprising all connected and 
topologically distinct diagrams constructed from Go and the Coulomb potential v. 
As an alternative, Hedin || derived a set of exact integral equations that define 
£ as a functional of G. This constitutes a closure relation, which, in combination 
with Dyson's equation, allows a self-consistent algebraic determination of the Green 
function. As intermediate quantities, Hedin's equations 



£(1,2) = i /G(l,3)W(l+,4)r(3,2;4)d(34), (5) 
W(l, 2) = w(l, 2) + / W(l, 3)P(3, 4)v(4, 2) d(34) , (6) 



P(l,2) = -i /'G(2,3)G(4,2)r(3 ! 4;l)d(34) ! (7) 
r(l,2;3) = <S(1, 2)5(1, 3) + | | G(4, 6)G(7, 5)r(6, 7; 3) d(4567) (8) 

employ the screened Coulomb interaction W, the polarizability P and the vertex 
function T. The notation (1 + ) indicates that a positive infinitesimal is added to the 
time variable. 

In principle, this set of equations could be solved by iteration from a suitable 
starting point, such as Go, until self-consistency is reached. However, the occur- 
rence of a functional derivative in Eq. (||) prevents an automatic numerical solution, 
because it changes the mathematical expression for the integrand in each loop. In 
practice this means that the starting point must be chosen so close to the expected 
solution that self-consistency is reached after a very small number of iterations, for 
which the functional derivative can be evaluated analytically. For solids the Hartree 
Green function seems a good enough starting point, and solving Hedin's equations 
with £ = then produces the so-called GW approximation |J 

E GW r(l,2) = iG(l,2)W(l+ 2) (9) 

after one iteration. The screened interaction enters in the time-dependent Hartree 
or random-phase approximation, in which the polarization propagator is given by 

PRP A (l,2) = -iG(l,2)G(2,l). (10) 

The GW approximation can be regarded as a generalization of the nonlocal Fock or 
exchange potential E x (l, 2) = iG(l, 2)v(l + , 2) with dynamically screened exchange. 
In addition to Pauli's principle it includes polarization effects and therefore describes 
correlation between electrons with parallel spin as well as electrons with opposite 
spin. 

To be consistent with the iterative solution of Hedin's equations, the GW ap- 
proximation should be evaluated with the Hartree Green function Go, although 
in practice the corresponding Kohn-Sham propagator is typically used [p|-^2]|. A 
comparison with the exact expression ([5]) identifies two types of omissions: lack 
of self-consistency and neglect of vertex corrections. The first group comprises all 
terms that allow the self-energy to be written in the mathematical form (||) and 
the polarizability in the form ( |l0| ) with suitably defined effective propagators. Ver- 
tex corrections, on the other hand, stem from an expansion of the vertex function 
T. In this case the topology of the resulting diagrams forbids a reduction to the 
simple forms above. Some self-energy diagrams are shown in Fig. [l| Arrows rep- 
resent noninteracting Green functions, the screened interaction is indicated by a 
wiggly line. The first diagram represents the GW approximation. The second is a 
self-consistency term, because it can be absorbed into the first by an appropriate 
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Figure 1: Diagrammatic expansion of the self-energy. Arrows represent noninter- 
acting Green functions, the screened interaction is indicated by a wiggly line. 

renormalization of the Green function. This is not possible for the third diagram, 
however, which therefore represents a vertex correction. 

Suitable indicators for the assessment of corrections beyond the GW approxi- 
mation are their influence on the band structure, the spectral function and the total 
energy. The band structure is obtained from 

E n k = £nk + ( < f 2nk|S(E„ k )|<Pnk) • (11) 

If the single-particle potential V s is chosen to include the exchange-correlation po- 
tential of density-functional theory, as is typically done in practical calculations, 
then the self-energy must be replaced by the difference E(i?„k) — V^ c . The quasipar- 
ticle energies feature as prominent peaks in the spectral function, which is propor- 
tional to the imaginary part of the Green function. However, the spectral function 
also contains other resonances related to many-body excitations. For instance, a 
quasiparticle excitation may be accompanied by a series of satellites that correspond 
to the creation or destruction of additional plasmons. Finally, the total energy is 
given by the Galitskii-Migdal formula [[lJJ 

E=- f duj I d 3 r lim[w + /i(r)]ImG(r,r';ct>) , (12) 

where the one-body Hamiltonian h contains the kinetic energy operator and the 
external potential. 



3 Self-consistency 

Although Hedin's equations define the self-energy as a functional of the Green func- 
tion, which is in turn self-consistently determined by the self-energy through Dyson's 
equations, most practical calculations simply evaluate the expression (||) using the 
Green function Go of an appropriate noninteracting system. This approach is con- 
sistent with Hedin's derivation of the GW approximation and generally yields quasi- 
particle band structures in good agreement with experiments. Nevertheless, from 
a fundamental point of view this procedure has several shortcomings. First, the 
use of Go in the self-energy introduces an ambiguity, because the result of a GW 
calculation then depends on the choice of the single-particle potential V s . Aulbur, 
Stadele and Gorling [jl5[ recently investigated this effect by comparing a GW cal- 
culation based on the standard local-density approximation with one starting from 
an exact exchange and local-density correlation potential. Although the substantial 
difference in the initial Kohn-Sham band gap is largely levelled out and reduced to 
about 0.1 eV for most materials considered, in some cases important quantitative 
discrepancies remain. For instance, for GaAs unreconciled band gaps of 1.16 eV and 
1.90 eV are obtained in this way, compared to the experimental value 1.52 eV. Sec- 
ond, without self-consistency the GW approximation fails to conserve the particle 
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number, energy and momentum under time-dependent external perturbations |lf 
Even in equilibrium the integral over the spectral function 
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N =-J du J d a rlmG(r,r;u) (13) 

does not equal the number of electrons JlTj , as it should, although the quantitative 
error is less than 1% for typical semiconductors |l8| and the homogeneous electron 
gas in the range of metallic densities |l9). Furthermore, different methods of calcu- 
lating the total energy from the Green function are not mutually consistent |2(J . 

Due to the high computational cost, fully self-consistent GW implementations, 
in which the Green function obtained from Dyson's equation (Q) is used to update 
the self-energy until convergence is reached, became possible only a few years ago. 
The first calculation was probably performed by de Groot, Bobbert and van Haerin- 
gen ]2jj for a quasi-one-dimensional semiconductor, in which the crystal lattice is 
modelled by a sinusoidal potential. The most important finding is a substantial 
increase in the band gap, which overestimates the exact Monte-Carlo result and 
comes close to the Hartree-Fock gap. Furthermore, the weight of the incoherent 
background in the spectral function is drastically reduced, leading to sharper and 
more pronounced quasiparticle resonances. 

Although the physical characteristics of the model raised doubts concerning 
its relevance for real materials [Q], recent calculations for silicon within a finite- 
temperature approach have confirmed these findings [^2| : while the standard non- 
self-consistent GW approximation widens the indirect band gap from 0.56 eV in the 
local-density approximation to 1.34 eV, in good agreement with the experimental 
value 1.17 eV, self-consistency increases the gap to 1.91 eV. The self-consistent GW 
result hence overestimates the experimental band gap by as much as the local- 
density approximation underestimates it. Besides, self-consistency again leads to 
an accumulation of spectral weight in the quasiparticle peaks, in disagreement with 
experiments. The quasiparticle peaks are also narrowed, corresponding to increased 
lifetimes. 

The majority of applications have so far focussed on the homogeneous electron 
gas, taking advantage of mathematical simplifications due to the spatial isotropy. In 
addition to fully self-consistent results (lj],^o|,^3],0, several studies have reported 
partially self-consistent calculations in which the Green function G in E q. (j9|) i s 
updated until convergence but the screened Coulomb interaction W is not ||25|-|27| . 
In contrast to conventional GW calculations, which give an accurate account of 
the correlation-induced band narrowing, self-consistency causes the occupied band 
width to increase even above its free-electron value fcp/2, where denotes the 
Fermi wave vector. Incidentially, an expansion of the valence band width was also 
found for silicon Q . The weight of the quasiparticle peaks is again increased and 
that of the plasmon satellite reduced accordingly. While the quasiparticle peaks 
are narrowed, the satellite is broadened and shifted towards the Fermi level. A 
calculation for potassium shows that the results of the homogeneous electron gas 
can be fully generalized to metals |2^] : the Kohn-Sham band width of 2.21 eV is 
narrowed to 2.04 eV in the first-order GW approximation, improving agreement 
with the experimental value 1.60 eV. The remaining discrepancy can actually be 
explained in terms of measurement effects, which shift the apparent peak position 
of the Fermi level more towards lower binding energies than the state at the bottom 
of the band and thus give rise to an additional artificial narrowing between 0.2 eV 
and 0.4 eV p8fl . In comparison, the self-consistent GW band width is 2.64 eV, much 
larger even than the Kohn-Sham value. 

Taken together, the above results form a consistent picture of the effects of self- 
consistency, which in a dramatic way reverses the correct trends of conventional 
GW band-structure calculations and destroys the originally good agreement with 
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experiments. This can be understood as follows. In contrast to Go, which refers to 
a noninteracting system with single-particle excitations only, part of the spectral 
weight in the interacting Green function G is transferred to plasmon satellites, which 
describe collective excitations. The quasiparticle peaks are reduced accordingly. 
This redistribution of spectral weight in turn implies a smaller dynamic self-energy, 
i.e., the part E c = iG(W — v) of the self-energy that is due to correlation, in the 
vicinity of the quasiparticle position. The dynamic self-energy is positive, tending 
to reduce the band width, and competes with the exchange part E x , which is always 
negative and increases the band width, e.g., by an amount k-p/ir for the homogeneous 
electron gas in the Hartree-Fock theory In a self-consistent calculation the 
smaller dynamic part no longer dominates over the exchange part, which is only 
marginally reduced, and the band width hence grows rather than narrows [p6| . The 
reduced dynamic self-energy also explains the larger renormalization factors 

<9Re£ c (fc,w) \ 1 
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(14) 



which indicate the weight of the quasiparticle resonances, as well as the increase 
in their lifetime, which is inversely proportional to the imaginary part of S c . The 
effects described above are further reinforced if the screened interaction is also 
calculated self-consistently [^0|. In this case the sharp plasmon excitations in W 
disappear and it no longer has a physical meaning as a response function. 

In conclusion, self-consistency is not a good idea for calculating quasiparticle 
energies without the simultaneous inclusion of vertex corrections. However, an 
interesting and surprising outcome of the self-consistent GW calculations for the 
homogeneous electron gas was the realization that the total energy obtained from 
the Galitskii-Migdal formula (|l2|) is strikingly close to supposedly exact Monte- 
Carlo data [Q. It has been speculated that this quite unexpected result is 
related to the fact that the self-consistent GW approximation is conserving in the 
sense of Baym and Kadanoff , although relative energy conservation under time- 
dependent external perturbations does not necessarily imply an accurate total en- 
ergy on an absolute scale. 

Based on the investigation of finite Hubbard clusters, Schindlmayr, Pollehn and 
Godby ^| note that self-consistency systematically raises the total energy due to 
an absolute shift of the chemical potential towards higher energies. Likewise, for 
the homogeneous electron gas the upward transfer of spectral weight from low-lying 
plasmon satellites to the quasiparticle peak and the increase in the band width, 
which moves the quasiparticles to lower energies relative to the chemical potential, 
work in different directions and largely cancel. On the other hand, the chemical 
potential is shifted upward on an absolute scale, explaining most of the change 
in the total energy. Although the increase in total energy often leads to better 
agreement with exact numerical solutions for the finite Hubbard clusters, in some 
parameter ranges the results become worse. In combination with the unphysical 
features of the spectral function described above, these observations suggest that the 
good quantitative agreement with Monte-Carlo data for the homogeneous electron 
gas may be fortutious. Recent calculations for the spin-polarized and the two- 
dimensional electron gas indeed show slightly larger errors [H. Nevertheless, these 
results have inspired renewed interest in total-energy calculations within many-body 
perturbation theory 



4 Vertex corrections 

Vertex corrections introduce additional interaction channels not accounted for in 
the GW approximation (^) for the self-energy and the random-phase approxima- 
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tion (|10|) for the polarizability. Their effect can be understood by physical inter- 
pretation. The random-phase approximation describes dynamic screening within a 
time-dependent Hartree approach. The screening electrons around a photoemission 
hole are thus considered independent, exchange and correlation, which enforce spa- 
tial separation, are ignored. As a result, the negative charge cloud is too tightly 
drawn around the central hole and screening at small distances is overestimated, 
so much, in fact, that the pair distribution function in the random- phase approxi- 
mation even becomes negative, which is unphysical ]32] ]. Vertex corrections in the 
polarizability will therefore, in general, reduce the screening and strengthen the in- 
teraction. They can also introduce new physical phenomena, such as bound states 
between the electrons and holes created during a photoemission process. In this 
case the response function acquires an additional exciton resonance that lies below 
the plasmon energy. Vertex corrections in the self-energy describe the same ex- 
change and correlation mechanisms between the central photoemission hole and the 
surrounding particles. The two effects compete and cancel partially: while vertex 
corrections in the self-energy reduce the probability of finding other holes near the 
central photoemission hole, vertex corrections in the polarizability simultaneously 
reduce the screening, thereby increasing the hole density. The change in the quasi- 
particle energies will hence be small, albeit important. On the other hand, the 
effect on the satellite spectrum may be drastic, because new types of excitations 
come into play. 

To be consistent, the same vertex function should be used in the self-energy and 
the polarizability |2^,^3|,^4). However, as the cost of calculating Feynman diagrams 
grows very rapidly with the topological complexity, nondiagrammatic vertices or 
plasmon-pole models for the screened interaction have long been the only way to 
determine higher self-energy terms. Ummels, Bobbert and van Haeringen [ pE[ eval- 
uated the lowest-order vertex correction displayed in Fig. |l| with a plasmon-pole 
model for silicon and diamond, revising an earlier calculation pq |. While the GW 
approximation increases the direct gap at the T-point by 0.78 eV for silicon and 2.12 
eV for diamond relative to the local-density approximation, in good agreement with 
experiments, the vertex diagram yields an additional contribution of —0.26 eV and 
—0.09 eV, respectively. The results for other high-symmetry points in the Brillouin 
zone are similar. It hence appears that the vertex corrections can be numerically 
significant, and the reduction in the band gap tends to cancel the increase due to 
self-consistency P2|| . However, the use of a plasmon-pole model leaves some uncer- 
tainties, and it is also known from an early application to the homogeneous electron 
gas that the lowest-order vertex correction leads to unphysical analytic properties 
in the self-energy, which can produce regions with a negative density of states j37| . 
These unphysical features are only cancelled by higher-order terms. 

When vertex corrections are taken into account both in the self-energy and 
the polarizability, their effects tend to cancel to a large degree. For the homo- 
geneous electron gas, Mahan and Sernelius |m| calculated the band width in the 
range of metallic densities with a static vertex function and obtained essentially 
the same result as in the GW approximation, although there is a strong addi- 
tional narrowing if vertex corrections are included only in the polarizability. The 
latter point was actually noted before in relation to the band structure of the al- 



kali metals [ll|l2|. Likewise, with a vertex function r(l,2;3) = 8(1, 2)/ xc (l, 3), 
where / X c(l,3) = 5V XC (1) / Sn(3) denotes the exchange-correlation kernel in the 
local-density approximation, Del Sole, Reining and Godby [Q found no change 
in the band gap of silicon, although the absolute position of the bands is shifted 
by about 0.4 eV. This vertex function is obtained by a consistent first-order solu- 
tion of Hedin's equations starting from density-functional theory Q as opposed to 
Hartree theory, which produces the GW approximation. Stronger changes occur 
if the vertex function is included only in the polarizability and not balanced by 
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corresponding self-energy diagrams: direct gaps are reduced by up to 0.2 eV and 
the valence band width decreases by 0.55 eV. Similar conclusions were reached for 
a finite Hubbard cluster |3!| . 

As expected from the generally good results obtained in the standard GW ap- 
proximation, vertex corrections, if applied in a consistent manner, tend to cancel, 
inducing subtle but quantitatively small changes in the band structure of solids. 
Many authors have also reported cancellation between combinations of vertex and 
self-consistency diagrams for a variety of systems [^5],|3f], This observation is 
reassuring, because it strengthens the theoretical foundation of the GW approxi- 
mation. On the other hand, it leaves the question unanswered how improvements 
might be achieved. Schindlmayr and Godby ||4lf proposed a systematic approach 
based on a continued iterative solution of Hedin's equations. As the GW approx- 
imation obtained after the first iteration represents a significant and systematic 
improvement over the zeroth-order Hartree or local-density-approximation, further 
advances might be achieved by a continuation of this procedure. The implicit inte- 
gral equation (||) for the vertex function can in fact be solved, and the functional 
derivative with respect to the full Green function G is at the same time replaced 
by a derivative with respect to the known propagator Go, which can, in princi- 
ple, be solved at all levels of iteration. The vertex correction obtained at the end 
of the second iteration starting from Hartree theory mixes diagrams of different 
order in the screened interaction. Numerical results for a finite Hubbard cluster 
show signs of convergence in the excitation energies and suggest that an iterative 
solution of Hedin's equations may improve the spectrum. However, this approach 
shares the problem of possible incorrect analytic properties that was earlier noted 
for the diagrammatic expansion of the self-energy by orders of the screened inter- 
action |37|. Takada j|2| proposed a similar approach in which Dyson's equation is 
solved self-consistently in each iteration, but so far it has only been exploited to 
derive a solution method for model systems without energy dispersion [ ff3[ . 

Incidentially, more progress has been made with respect to systematic improve- 
ments of the satellite spectrum, which from the start is not well rendered in the GW 
approximation. In particular, for the homogeneous electron gas and the alkali met- 
als, the GW approximation only produces a single plasmon resonance instead of a 
sequence of satellites separated from the main peak by multiples of the plasmon fre- 
quency. The so-called cumulant expansion remedies this problem by describing the 
coupling of a quasiparticle to multiple plasmons |44| . It includes the vertex diagram 
displayed in Fig. [l as well as corresponding higher-order terms. An application to 
Na and Al significantly improves the agreement with experimental photoemission 
spectra, although the relative intensities of the satellites with respect to the quasi- 
particle peak are still in discrepancy f45f| . The quasiparticle position is the same 
as in the GW approximation. In this context it is interesting to note that with 
vertex corrections, the self-energy is less affected by the deterioration of spectral 
features, such as the accumulation of spectral weight in the quasiparticle peak and 
the broadening of plasmon satellites, when applied self-consistently ]27|] . This is 
another indicator that the vertex corrections contain the correct physical features. 
If the exact vertex function is used, then self-consistency must, of course, yield the 
true spectrum without the mentioned adverse effects. 

In a similar way, excitonic effects in the optical absorption spectrum of semicon- 
ductors can be accounted for by vertex corrections that describe multiple scattering 
or binding between electron-hole pairs. This procedure is very computationally 
demanding, because it requires a calculation of the two-particle Green function. 
Nevertheless, excitonic vertex corrections have been calculated for several semicon- 
ductors and significantly improve the agreement with experimental spectra 
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5 Summary 



The GW approximation is a reliable method for ab initio electronic-structure cal- 
culations that produces band structure in good agreement with experiments for a 
wide range of materials. This was originally taken as a sign that higher-order self- 
energy diagrams are negligible. Statements to this effect indeed frequently appear 
in the early literature Q] . Only the availability of modern computational resources 
in recent years has it made possible to explicitly evaluate self-consistency and vertex 
corrections and test this assumption. The numerical results show that the excluded 
terms are, in general, not individually small but tend to mutually cancel. The find- 
ings of self-consistent calculations are consistent and show a serious deterioration 
of spectral features compared to the standard GW approximation, which can be 
understood by shifts of spectral weight in the dynamic self-energy. The influence of 
vertex corrections is naturally less clear-cut, reflecting the large variety of possible 
vertex functions. However, the case for a mutual cancellation of vertex corrections 
in the self-energy and the polarizability is well established. Although a general and 
systematic way of improving quasiparticle energies is still outstanding, physically 
motivated vertex corrections for better satellite spectra are known. 
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